function [quad_w, quad_xy] = quad_rule_dunavant(d)
% function [quad_w, quad_xy] = quad_rule_dunavant(d)
% 

if (d > 20 || d < 1)
    error('Sorry, I can"t support the degree large than 20!\n');
end
fname = strcat('dunavant/dunavant_quad_tri_',num2str(d));
load(fname,'xy','w');
quad_w = w; % the original is for large triangle
quad_xy = xy;